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Сформульовано обернену контактну задачу гравіметрії для контактної поверані в класі 
Сретенського. Доведено розривність вертикальної пожлідної сили тяогюіння, яка фігурує 
в її правій частині. Цю особливість було використано для редукції задачі до альтер- 
нативної постановки у вигляді інтегрального рівняння 1-го роду. Шлятом різницевого 
аналізу меою відповідних інтегралів виведено ефективний спосіб послідовного обчислення 
контактів. Задача узагальнюється на п. меос та зберігає коректність у межах класу 
Сретенського. 


З початку ХХІ ст. відбувається неперервна зміна основної парадигми теорії інтерпретації 
потенціальних полів ||. Вона супроводжується переглядом ключових методологічних |З, 
апаратно-програмних |З| та числових |4| засад тлумачення даних цих полів. У рамках ство- 
рення числового базису нової парадигми особливого значення набуває розробка адекватного 
геофізичній практиці апарату математичного моделювання геофізичних подів, адаптовано- 
го до обробки масивів даних великої розмірності, в тому числі й на комплексах паралельних 
обчислень. 

Ключове місце в рамках нової методології належить аналітичним апроксимаціям се- 
редовища й поля. Теоретичне підгрунтя такого підходу викладено в праці Страхова |З. 
У цьому руслі активно розробляються конструктивні аналітичні апроксимації геологічного 
середовища. Зокрема, введені в науковий обіг апроксимації на основі типових блоково-шару- 
ватих або сферично-циліндричних елементів |5|. Постульовано нагальну потребу доповнити 
ці конструкції апроксимаціями потенціальних полів |б|. Це, певною мірою, було зроблено 
в статті |7|, де запропоновано новий математичний базис (диференціальне рівняння) для 
опису аномалій сили тяжіння. Також вивчаються апроксимації рельєфу поверхні спостере- 
жень шляхом параметризації заданими гладкими функціями |5Ї. 

Продовжуючи актуальні розробки аналітичних апроксимацій середовища, ми пропонує- 
мо нову апроксимаційну схему для визначення контакту двох однорідних тяжіючих тіл, 
які належать до класу Сретенського, за заданим розподілом вертикальної компоненти ано- 
малій сили тяжіння. Зауважимо, що аналітичні конструкції методу дозволяють реалізувати 
розпаралелювання обчислень за рахунок рекурентності виразів для розв'язання відповід- 
ної прямої задачі. 

Відзначимо, що постановка та розв'язання прямих 1 обернених задач граві- й магніто- 
метрії у фомулюванні для класу тіл Сретенського започатковані в працях Є.Г. Булаха |З. 
Ці задачі викладено для числової моделі як з одного |10Ї, так і кількох |11| рудних тіл, У си- 
лу вдалої та зручної параметризації вони зручні для практичного застосування в комплексі 
інтерпретації даних рудної геофізики. Однак формулювання самого класу Сретенського 
моделей аномальних тіл дозволяє представити пряму задачу гравіметрії у структурній пос- 
тановці, чого дотепер не було зроблено. Поширення розв'язку цієї задачі на кілька кон- 
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Рис. 1. Модель рудного тіла, для якої середня площина у смислі Л. Сретенського є неоднозначною 


тактних поверхонь відповідатиме сучасній практиці моделювання складного геологічного 
середовища за комплексом геофізичних даних. 

Постановка задачі. Нехай на деякій горизонтальній площині 2 - 0 (вісь 2 спрямована 
вниз) задані значення вертикальної похідної Г, (7, 0) потенціалу сили тяжіння, зумовленого 
"нескінченним" шаром С (подовжується в обидва боки за профіль спостережень), видовже- 
ним уздовж осі у. Нехай шар С складається з двох пластів із постійною густиною с з соп5і, 
розділених горизонтальною контактною поверхнею Н. Цей об'єкт належить класу двови- 
мірних тіл, для яких існує середня площина 20 (вона проходить через тіло таким чином, 
що будь-який перпендикуляр до неї перетинає поверхню тіла лише в двох точках, по різ- 
ні сторони від площини). Контактна поверхня Н ундулює навколо середньої площини 70, 
не надто ухиляючись від неї, у смислі метрики банахового простору | Н - Д0ЇЇв (? р 


У частинному випадку Н збігатиметься із 29. Назвемо цей модельний клас тіл -- класом 
Сретенського Зг(1, Д), (структурний варіант), де 1 - - означає сталу густину с, а 0 - замк- 
нену область (поверхню), виповнену тяжіючими масами. 

Зауваження 1. Якщо площина спостережень (денна поверхня) не є горизонтальною, Її 
можна, апроксимувати будь-якою конструкцією 7 - Н(2, у), в тому числі однією із пропо- 
нованих в статті |З. 

Зауваження 2. Визначення класу Зг(1, Д) через поняття середня площина не є конструк- 
тивним, оскільки застерігає лише від вжитку "глечикоподібних" тіл і ніяк не регламентує 
спосіб визначення середньої площини, наприклад, для множини багатокутників (рис. 1). 

На нашу думку, подолати амбівалентність означення класу Зг(1, 0) можна, залучивши 
диференціальні характеристики конкретного сімейства тіл, як, наприклад, зроблено в |З. 
Утім, подібні міркування виходять за рамки даного повідомлення. Модель горизонтального 
циліндра апріорі задовольняє вимоги класу 8г(1, Р). 

Уточнимо характеристику класу 8г(1, 0). Поєднуючи із середньою площиною аномаль- 
ного тіла координатну площину со), аналітично опишемо область Р, зайняту тілом 


ра(8О:а«є«БСЯТЄ) є с « СО), (1) 


де -со «С (0) « осо -- певні контактні поверхні. При цьому межу 00 області 0) можна 
представити як об'єднання двох контурів 01 | | 012»: 


92р; - (5бО:а«є«вСс- сі - 1,2). (2) 


Розглянемо підклас 5г(1, До) класу 8г(1, Р), коли а - - 00, - о0, аосі координат Є й 5, С 
й 2 паралельні. Тоді, не порушуючи загальності міркувань, можна представити область Що 
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та її межі в аналітичному вигляді: 
роз (8,0: ооо ою, 2) Є С є 29), 
9Дої З ЦО: - оо «Є « со, 20) «з СОС), і- 1,2). 
Нехай підклас 5г(1, До) тяжіючих тіл складається з двох контактних поверхонь 
І «сф, | -со«сас, | із-1,2, (4) 


розділених середньою площиною 2 з 2) » 0. 

Поставимо задачу: визначити контакти (4) за значеннями вертикальної похідної Г, (2, 0) 
потенціалу сили тяжіння. 

Аналітичні властивості аномалій. Розв'язати вказану задачу допоможе аналіз влас- 
тивостей похідних потенціалу сили тяжіння Ї/(х, 2). Представимо ці похідні як 


аа) -9ге Шен ор 
(5) 
ПОС. ль 1 су 
де 
ро: оче « соб «єс (6) 


Внаслідок подання формули (6) таким чином, подвійні інтеграли в формулах (5) можна 
переписати у вигляді 


с) 
4- 
пам їі «є і 
0 
с) 


пн Ї є пед 


Зважаючи на результати, що отримані в статті |12|, після ряду нескладних аналітичних 
перетворень матимемо співвідношення 


п КТ 
9С (2, 2) т (є - а)2-6 СО) 2 сС (2), з. І); (7) 
а зару -а 

дз -г| ер оо) ЗміеїКо) с да, Ока, 


-2п оС (а), (т) є 2. 


Отже, функція М, (2,2) щодо 2 є неперервною, в той час як функція (2,2) щодо 2 
є розривною. На підставі дослідження фундаментальних властивостей (7) контактів (4), 
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згідно з методикою аналізу гармонічних функцій, розвиненою в публікації |12|, отримуємо 
при 2 -» - 0 такі вирази зовнішніх похідних гравітаційного потенціалу: 


2 


дл) 2. -мУр | мав зАоиМ 


2 


да 2) їі (Є -- аг)? (С (ау 
рога ше гУіа | іс - | й 





При цьому на середній площині 2 -- 29 (С9)(2) « 2 « С(а), -о0 « ж « оо) з тих самих 
співвідношень (7) знаходимо, що 


9С (2, 20) - ї СО - а 
У о 3. іі агсіє ої 


2 БУ 
9С (2, 2) 1 2 (2 -а)2-н|С9У(2)|? 
ОО дпа С) 2)-2л Й» сФ 2)-220|Ї СЕ а 
92 (2) С" (2) -220| Ж ве 
1-1 За 

Спосіб розв'язання. Очевидно, отримано розбіжності значення потенціалу на поверх- 
ні (8) і на середній площині (9). Це нескладно було передбачити, виходячи з співвідношен- 
ня (7). Однак у цьому протиріччі закладена цікава можливість для розв'язання поставленої 
задачі. 

Дійсно, знайдемо внутрішню межу похідної при обчисленні інтеграла Пуассона: 
ги. їду 

(2,0 1 уз 2 

М ) рАро і (п о) 0 54 (10) 

да т ді  (п-а)2- 2 


-оор 


. Представимо підінтегральну функцію у вигляді 


сно) 2). а іт Пе шу заейС За 1 о 


для тяжіючої області БД; с 1(6,С): - 00 « Є « ос0,0 «(С « сф (2)). Обчислюючи інте- 
грал (10), отримуємо з урахуванням наведеного вище зображення такий вираз: 





"б 2 т. ад дій 
п з панна а) СЛ ї што 11) 


Враховуючи значення зовнішньої межі (8) та виразу (11), знайдемо різницю меж 


да 0) да, 0 Ї. (є-а)?а КФ(дур 
че Е го нама МАрУнн 1 її - -а)зя ст М 
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Щоб спростити подальші викладки, позначимо цю різницю потенціалів через функцію 


М (а): 


Відтак, нами отримано нелінійне інтегральне рівняння 





17. (Є - а)? (СО 
(Фа) з / пото М 12 
ія ро Ма 
для визначення першого наближення контакту 20) С (г), -ор «ж « оо. Розв'язання 


подібних рівнянь започатковане в публікації |12|. 
Визначивши із останнього рівняння функцію (1) (2), можна обчислити наближення поля 


оо 


- -дтренсі а) - ор (Пи 


-оо 


оС (ов 0) 
02 


Кей зн ар 


Є - а)? ФО г 





А це, в свою чергу, дозволяє знайти в "чистому вигляді" ефект від подальшого контак- 
ту -- у вигляді різниці граничного та попереднього значень вертикальної похідної потенціалу 
сили тяжіння: 

905 (х,0)  д0,.(х,0) дО (2,0) 
02 02 дх 
Маючи цю величину, можна з рівняння, аналогічного рівнянню (13), тобто 


оо 





- є) |СЧа)Ї 005 (ж, 0) 
заблніуно ас о ЗР) а 
обчислити чергове наближення контакту 22) С (2) (2), -00 « г « оо. Таким чином, було 


знайдено "верхню" і "нижню"? межі області (1), що й розв'язує поставлену задачу (4). До 
речі, власне, самі контакти С (9 (Є) можна апроксимувати конструкціями, описаними в |11Ї. 

Зауваження 3. Наведені вище математичні викладки легко узагальнюються на випа- 
док п меж 20) -- С?) (2), -00 «0 22 «0 00, якщо між кожною парою цих меж можна провести 
середню площину (тобто якщо ці межі не виходять за межі класу 8Г(1, Д)). 

Таким чином, нами введено в науковий обіг нові апроксимаційні конструкції (12) й (13), 
(14), за допомогою яких можна ефективно відновити будову шаруватого геологічного се- 
редовища, зображеного парами 2 -- С?) (2), і з- 1,2 (-верхня" і "нижня"? межі) контактів, 
обмежених властивостями класу Сретенського 5г(1, Р). 

Доведено, що на класі 5г(1, Р) послідовні наближення С? (г) однозначні й стійкі (збі- 
гаються до заданого точного розв'язку (10) (х)), якщо похибка вхідних даних не переви- 
щує 390 від максимальної амплітуди аномалій М, (т, 0) |13|. 

Дискретизацію інтегральних рівнянь (12) й (13), (14) доцільно здійснити за допомогою 
відомих методів скінченновимірної редукції інтегральних рівнянь 1-го роду |14|. Алгоритм 
розв'язання задачі пройшов первинну апробацію на працездатність. У перспективі передба- 
чається поширення аналітичних викладок на тривимірний випадок та апробація алгоритму 
на практичних матеріалах. 


155 1025-6415 | Доповіді Національної академії наук України, 2015, Мо 3 91 


10. 


11. 


12. 


13. 


14. 


Стралов В. Н. Принципиально новая теория интерпретации данньх о потенциальньтх полях (грави- 
тационньх и магнитньх аномалий) // Геофиз. журн. - 2003. - 25, Хе 1. - С. 3-7. 


Стразлов В. Н. Об зффективньгх по бьістродействию и точности методах построения линейньгх ана- 
литических аппроксимаций в геофизике, геоинформатике и гравиметрий // Там же. - 2007. - 29, 


Хе 1.- С. 56-84. 


Бьшчков С. Г., Симанов А. А., Хохлова В. В. Современнье процедурві вьічисления аномалий сильт 
тяжести при вьгсокоточньх гравиметрических наблюдениях / / Вестн. Перм. ун-та. Геология. - 2013. - 


Вьш. 3(20). - С. 61-70. Беєр:/ /суБегіепіпКа. ги /агбісіе/п /. 


Стралов В. Н., Керимов Й. А., Степанова Й. З. и др. Новьшй информационньй базис гравиметриий и 
магнитометрим / / Геофизика и математика, - Пермь: Горн. ин-т УрРО РАН, 2001. - С. 274-277. 


Страхов В. Н Линейнье аналитическиєе аппроксимации рельефа поверхности Земли / / Геофизика и 
математика. - Москва: Изд-во Моск. ун-та, 1999. - С. 12-19. 


Якимчик А. ІЇ., Чорна О. А. Про побудову аналітичних апроксимацій елементів аномальних гравіта- 
ційних полів // Вісн. Київ. нац. ун-ту ім. Тараса Шевченка. Геологія. - 2010. - Вип. 51. - С. 12-14. 


Дубовенко Ю. І. Відновлення потенціалу сили тяжіння за значеннями модуля його градієнта в задачі 


Алексідзе // Там само. - 2011. - Вип. 55. - С. 61-65. 


Дубовенко Ю. И., Черная О. А. Об определениий плотностного контакта в сложнопостроенной среде / / 
ХІ Междунар. конф. по геоинформатике - теоретические и прикладньге аспектьгп ("Геоинформати- 
ка-20127), Киев, 14-17 мая 2012 г.: Сб. науч. тр. - Киев: ТОВ "Карбон ЛТД?, 2012. - СР-ВОМ. - 
С. 3494. 


Булах Е. Г. Обратнье задачи гравиметриий для геологических моделей класса Сретенского // Доп. 
НАН України. - 2002. - Хе 1. - С. 117-119. 


Булах Е. Г. О другом аппроксимационном построениий геологической модели класса Сретенского для 
решения обратньх задач гравиметрим / / Там само. - 2002. - Хо 5, - С. 128-132. 








Булах Е. Г., Слободник Н. А. Обратньге задачи магнитометрии для совокупности тел класса Л. Н. Сре- 
тенского / / Геофиз. журн. - 2008. - 30, Хе 3. - С. 49-55. 


Дубовенко Ю. І., Чорний А. В. Дослідження оберненої задачі потенціалу для контактної поверхні // 
Там же. - 2002. - 24, Хе 3.- С. 77-92. 


Дубовенко Ю. И. Об определениий плотностньх неоднородностей в классе Сретенского // ХУ Урал. 
молод, науч. шк. по геофизике, 24-29 марта 2014 г., Екатеринбург: Сб. докл. - Екатеринбург: ИГф 
УРрО РАН, 2014. - С. 84-86. 


Верлань А. Ф, Сизиков В. С. Интегральнье уравнения: методь, алгоритмьв, программь. - Киев: 
Наук. думка, 1986. - 544 с. 


Інститут геофізики їм. С. І. Субботіна Надійшло до редакції 07.10.2014 
НАН України, Київ 


92 


Ю. И. Дубовенко 


Об определений контакта в модельном классе Сретенского 


Сформулирована обратная контактная задача гравиметрий для контактной повераности 
в классе Сретенского. Доказана разрьвность вертикальної производной, силь, тягжюести, 
фигурирующей в єе правой части. Зта особенность ипспользована для сведения задачи к аль- 
тернативной постановке в виде шитегрального уравнения 1-го рода. Путем разностного 
анализа границ соответствующих шитегралов введен зффективньгй способ последователь- 
ного ввічисления контактов. Задача обобщаєтся на тп границ и сограняєт, корректность 
в рамка класса Сретенского. 
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Уп. І. ДрибоуепКо 


Оп бе аейпійоп ої а сопбасі утїбріп їре ЗгебепякКії тоає! сіаз55 


ТРе іпиуєтзе сотіасі ртобієт, о) дгачібу іпуетяїот Їотг а депзіду іпістІасе шійут Ше бтеветзкії сіаз8 із 
розей. А аізсопбітпийу о) пе утауйу фетбвіса! аетітабіуєе, шРісі, арреатз от ртобіеті з гідпі-Пата зіаєе, 
із ртгочеа.. ТРіз Теаїите 15 шзеа їо тейисе Пе ртобієт іо їпе айетпайіе зїаїететі іп їйе Їогт. ої Не 
17 ріпа іпісута! едиабіот. Ву театз о) їе тезійиа! апаїузіз о) Ре боипаатіся Рог їе соттезропаїту 
іпісутаіз, ап ефбістепі кесітідие Їот те зиссеззїце саїсиіабіот о пе астпзібу ітіетіасе 55 дечеЇорей. 
ТВе ртобієт, із депетайіхеа отіо Пе сазе ої п боцпаатієз, шііїе ртезегиїпд ії5 соттесітез8 шійніт, 
їйе бгеїепзкії сіаз58. 
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